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Abstract 



We study electrically charged asymptotically flat black brane solutions whose world- volume fields 
are slowly varying with the coordinates. Using familiar techniques, we compute the transport 
coefficients of the fluid dynamic derivative expansion to first order. We show how the shear and 
bulk viscosities are modified in the presence of electric charge and we compute the charge diffusion 
constant which is not present for the neutral black p-brane. We compute the first order dispersion 
relations of the effective fiuid. For small values of the charge the speed of sound is found to be 
imaginary and the brane is thus Gregory-Laflamme unstable as expected. For sufficiently large 
values of the charge, the speed of sound becomes real, however, in this regime the sound mode 
attenuation is found to be negative. The electrically charged brane is thus found to be (classically) 
unstable for all values of the charge density in agreement with general thermodynamic arguments. 
Finally, we show that the shear viscosity to entropy bound is saturated, as expected, while the 
proposed bounds for the bulk viscosity to entropy can be violated in certain regimes of the charge 
of the brane. 
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1 Introduction 



Black branes possess hydrodynamic properties in addition to their thermodynamic properties. 
Indeed, fluid dynamics is the natural local generalization of the global thermodynamics of black 
branes. This is well-known and established in the hydrodynamic limit of the AdS/CFT corre- 
spondence [1], but has also been shown to hold true for Ricci flat black branes [2, 3, 4, 5]. Here the 
effective dynamics of the black brane (collectively known as 'the blackfold approach') is captured 
by the quasi- local stress tensor of Brown and York [6]. In general, the stress tensor is given by 
a fluid-elastic derivative expansion in the collective parameters describing the brane. To leading 
order, the stress tensor is that of a boosted perfect fluid which accounts for the thermodynam- 
ics of the brane. Fluctuations in the directions transverse to the brane (i.e. bending) give rise 
to elastic contributions [7, 8, 9] while fluctuations in the world-volume fields longitudinal to the 
brane, relevant to this paper, are captured by fluid dynamic dissipative corrections to the effective 
stress tensor [4]. 

In this paper we examine how the presence of electric Maxwell charge modifies the fluid 
transport coefficients of the neutral black brane originally considered in [4]. Moreover, we compute 
the transport coefficient associated to charge diffusion in the charged black brane to first order in 
the fluid derivative expansion. Although many supergravity solutions carrying lower dimensional 
charge on their world-volume are known [10, 11], only recently a new family of black branes of 
Einstein gravity coupled to a dilaton and a single q+1 form gauge field was found [12]. This family 
of solutions has both the dilaton coupling and the (integer) dimension q as free parameters. The 
Maxwell charged black brane solution can therefore straightforwardly be obtained from this family 
of solutions by specializing to q = and turning off the dilaton coupling. This especially means 
that we can study hydrodynamic fiuctuations of pure Einstein gravity coupled to a single gauge 
field without a dilaton - something which is not possible for generic supergravity backgrounds. 
In this way, our computation provides the simplest generalization of the neutral case. In detail, 
we consider long wavelength fluctuations around the black brane solution of Einstein-Maxwell 
gravity (which we shall from now on dub the Reissner-Nordstrom black brane) following the 
method of [1, 4]. We solve the full set of coupled Einstein-Maxwell equations to flrst order in the 
derivative expansion and compute the effective stress tensor and current. This provides us with 
the charged generalizations of respectively the shear and bulk viscosity along with the charge 
diffusion constant which is not present in the neutral case. 

Having computed the shear and bulk viscosities we flnd that the bound rj/s > l/47r, originally 
proposed in [13], is saturated. This agrees with the expectation that this should hold for any 
two-derivative gravity theory [14, 15]. However, we flnd that the bulk viscosity to entropy bound 
proposed in [16] can be violated in certain regimes. Perhaps this is not too surprising since the 
derivation of this bound relies heavily on holographic considerations. However, it is worth noting 
that the bound is saturated for the neutral black brane [4]. Finally we mention that the modified 
bound proposed in [17, 18] is found to be violated. 

By computing the speed of sound in the effective fluid of the neutral black brane ref. [3] was 
able to identify the unstable sound mode of the effective fluid with the Gregory-Laflamme (GL) 
instability [19, 20]. This (very simple) computation can already be carried out at the perfect 
fluid level. The results of [4] allowed further refinement of this result and showed remarkable 
agreement with numerical data. Performing a similar computation for the Reissner-Nordstrom 
black brane, we find that for a given temperature, the speed of sound becomes real for above a 
certain threshold value of the charge density. This seems to imply that the Reissner-Nordstrom 
black brane is not GL unstable for sufficiently large values of the charge. However, upon closer in- 
spection of (next-to-leading order) dispersion relations, we find that above the threshold value of 
the charge the attenuation of the sound mode becomes negative. This means that, in this regime, 
sound waves are not damped but instead growing exponentially. The Reissner-Nordstrom black 
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brane is therefore leading order and next to leading order GL unstable below and above the 
charge threshold, respectively. This complementary behavior of the instability is expected from 
the general thermodynamics. Indeed, by computing the specific heat C and the isothermal elec- 
tric permittivity c we find that the conditions C > and c > put complementary conditions 
on the charge density of the brane (recall that C < for the neutral brane). Since thermody- 
namic stability requires that C > and c > [21], the Reissner-Nordstrom black brane is never 
thermodynamically stable. 

In many ways, studying intrinsic fluctuations of branes in the blackfold formalism is similar in 
spirit to the well-known fluid/gravity correspondence of AdS/CFT [1]. We also mention that the 
computation in fluid /gravity analogous to ours (fluctuations of the AdS Reissner-Nordstrom brane 
of co-dimension 1 with a Chern-Simons term) was carried out in the papers [22, 23]. However, 
we emphasize that our computation deals with asymptotically flat branes of general co-dimension 
(note that extensions of the blackfold formalism to AdS backgrounds have been made in certain 
regimes - see [24, 25, 26, 27]) and that the effective fluid stress tensor has no direct interpretation 
as a "dual" fluid of a QFT. Also note that the blackfold fluid stress tensor is not that of a 
conformal fluid. Recently a connection between the fluid/gravity correspondence and the blackfold 
formulation was established. This was done by constructing a map from asymptotically AdS 
solutions compactified on a torus to a corresponding Ricci-flat solution obtained by replacing the 
torus by a sphere [28] . This was used to take the general second order results of fluid / gravity [29] 
and map them to the second order blackfold stress tensor. This provides even further improvement 
of the dispersion relation of the GL instability. 

The outline of the paper is as follows: In the subsequent section 2 we will setup the leading 
order solution and review its thermodynamics. In section 3 we discuss the perturbation procedure 
and explain how the boundary conditions are handled. In section 4 the flrst order equations are 
solved and in section 5 the effective stress tensor and current are provided. In section 6 the 
transport coefficients are used to analyze the dispersion relations and the GL instability. Finally, 
we will discuss the results in section 7 and address some interesting future developments. 

Notation: We use fi, v to label the D = p -|- n -|- 3 spacetime directions. Moreover, we denote 
the p + 1 world-volume directions of the brane in Schwarzschild coordinates by = (t,x*) and 
in Eddington-Finkelstein coordinates by a" = (u,(T*) with a = 0...p and i = The 
co-dimension of the brane is n + 2. For simplicity of the presentation we restrict ourselves to the 
cases where n > 1 due to a slightly different behavior at infinity for the n = 1 solution. However, 
treating the special case of n = 1 should be straightforward using similar considerations as for 
the neutral case. 

2 Reissner-Nordstrom branes and effective zeroth order fluid 

In this section we review the generalized Gibbons-Maeda solution for q = which was found in 
[12]. The generalized Gibbons-Maeda solution describes a black p-brane with horizon topology 
S"-~^^ X W which has electric g-charge diluted on its world-volume. The solution was obtained 
from the Gibbons-Maeda solution [30] through an elaborate double uplifting procedure. The 
general solution is given in terms of a metric, a dilaton and a (q + l)-form gauge field under 
which (the g-charge diluted on the) black p-brane is charged. A particularly nice property of 
the generalized Gibbons-Maeda solution is that the dilaton coupling a can be treated as a free 
parameter. This especially means that we are free to set a = 0. This is of course not possible 
for the well-known supergravity solutions such as the DO- Dp system [10, 11]. Moreover, in this 
paper we restrict ourselves to the q = case (Maxwell charge). 
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2.1 Reissner-Nordstrom black branes 

As explained above, we consider branes of Einstein-Maxwell theory. The action is 



1 



R - ^F^uF^'' 



(2.1) 



where Ffj_u is the field strength of the Maxwell gauge field A^, F = dA. We now present the 
Reissner-Nordstrom black brane solution. The solution is characterized by p flat spatial directions 
x*, a time direction t, and finally a radial direction r along with the usual transverse sphere S"""*"^. 
The total spacetime dimension D is related to p and nhjD=p + n + 3. The metric is given by 



\ i=l / 

The two harmonic functions / = /(r) and h = h{r) are given by-"^ 

f{r) = l-{^y, Mr) = l+(^)"7o. (2.3) 

The two parameters ro and 70 are related to the thermal and electrostatic energy of the solution 
(see below). The parameter B = B{p,n) is given by^ 

B = . (2.4) 

n+p 

Finally the gauge field A is given by 

A = V7o(7o + l)dt , (2.5) 

where we have defined N = B + 2. 

It is straightforward to apply a uniform boost to the solution (2.2) in the flat directions: 
The metric of the non-dilatonic boosted Reissner-Nordstrom black brane is given by 



ds^ = h"" \^-h-'^ f UaUbdx^dx" + f-'dr^ + r^dnf^^^^ + Aaidx^dx^j . (2.6) 

where A"^ = 5^^ + u°'Ub is the usual orthogonal projector defined by the boost u"". The gauge field 
is given by 

^ = ^ (7)" v'7o(7o + l) n^dx'^ . (2.7) 

In the next section we review the thermodynamics and effective blackfold description of this 
solution. 



2.2 Thermodynamics and effective blackfold fluid 

The blackfold theory of p-branes supporting q = charge was developed in [12] and further 
examined in [31] where also p = q was considered in detail (in a supergravity setting). For a 
uniform boost u"" of the brane, it is instructive to write the effective blackfold stress tensor in the 
form 

T^b ^ La^b _ }_^ab\ ^ ^Q^a^b ^ (2.8) 



^In the blackfold literature 70 = sinh^ a. 

^The full generalized Gibbons-Maeda solution has an additional parameter A. However, for the non-dilatonic 
Reissner-Nordstrom solution one has A — 2. 
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where ^ab is the induced metric on the blackfold. For our purposes (flat extrinsic geometry), 
we have 7^6 = r/afe- Moreover T is the local temperature, s is the entropy density, Q is the 
charge density and finally <1> is the electric potential conjugate to Q. The various quantities are 
parameterized in terms of a charge parameter 70 and the horizon thickness tq: 

The stress tensor in the form (2.8) immediately allows us to identify the thermal and the electro- 
static parts. Since Tq ~ Ts, tq gives us a measure of the thermal energy (density) of the given 
solution. In a similar manner 70 is identified with the thermodynamic ratio, 

1 $Q 

and 7o therefore measures the electrostatic energy relative to the thermal energy of the black 
brane. 

Of course, it is straightforward to cast the stress tensor into standard form (using the Gibbs- 
Duhem relation w = e + P = Ts + ^Q) 



T^-^ = wu^u^ + Pr?'^^ = eu^u^ + PA"^ , (2.11) 



where 



^-^^«"(" + ' + "~-)- ^-^^S. ^-=^^0"(l + A'70). (2.12) 
Finally the q = Q current supported by the p-brane is given by 

J(0) = Q?^"- (2.13) 

To leading order, the intrinsic blackfold equations take the form of the world- volume conservation 
equations SlaT^Q-^ = and Va J^q^ = 0. For flat extrinsic geometry 7^6 = -qab, they evaluate to the 
equations 

e = -w§^ = -w-^A'^'^dbP, Q = -Q^, (2.14) 

where = dau"" is the expansion of u'^ and a dot denotes the directional derivative along u". The 
(first order) equations will be important in the perturbative analysis. As expected, they will show 
up as constraint equations when solving the Einstein-Maxwell system perturbatively. 



3 The perturbative expansion 

As explained in the introduction, the aim of this paper is to solve the Einstein-Maxwell system 
in a derivative expansion around the solution given in section 2. In this section we define the 
appropriate coordinates to handle this problem and explain how the perturbations are classified 
according to their transformation properties under SO(p). 
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3.1 Setting up the perturbation 

Before perturbing the brane, we first need to cast the metric (2.6) into Eddington-Finkelstein-hke 
(EF) form. The reason is two-fold. First, it is essential for the computation that we can ensure 
regularity at the horizon and since the Schwarzschild description breaks down at the horizon, it 
is clearly more useful to use EF coordinates. Secondly, since a gravitational disturbance moves 
along null-lines, in order to control the perturbation, we want the lines of constant world-volume 
coordinates to be radial null-curves i.e. grr = 0. This is exactly the defining property of EF 
coordinates. For a general boost u", we define the EF coordinates cr" by 

a'' = x'' + u''n, r^{r)=r + j \ dr . (3.1) 

Here is chosen such that = h^^'^ / f and — )• r for large r. The first condition ensures 
that grr = while the latter is chosen such that the EF coordinates reduce to ordinary radial 
Schwarzschild light cone coordinates for large r. Notice that it is possible to write down a closed 
form expression for in terms of the hypergeometric Appell function Fi 



where the last equality applies for large r and is valid up to O ( ^.ai_i ). It is nice to note that the 
hypergeometric Appell function F\ reduces to the ordinary hypergeometric function in the 
neutral limit 70 — ?• 0. Indeed 

hm r,(r) = r,(r) = r + / Z':^') dr = r 2^1 f 1; -i; 1 - -; , (3-3) 

70^0 70 =0 Jr \ f J \ n n r^^J 

which is the used in [4]. With this definition of we will limit our analysis to the case for 
which n > 2. In EF coordinates, the metric (2.6) takes the form 

dsjo) = [-h~^ f UaUb da^da^ - 2h^^l^ Ua da'^dr + Aab da^da^ + r^dJ^J^+i)) . (3.4) 

Here the subscript indicates that the metric solves the Einstein-Maxwell equations to zeroth order 
in the derivatives. Notice that in these coordinates the gauge field will acquire a non-zero Ar 
component . However, we shall work in a gauge where this component is zero. We therefore take 



^(0) ^ /!;o \ " ^^0(^0 + l)'iiado-" , and in particular = , (3.5) 



h \ r 

Having determined the EF form of the metric and gauge field, we are now ready to set up the 
perturbative expansion. 

Following the lines of [4], we promote the parameters u"',rQ and 70 to slowly varying world- 
volume fields: 

u-^u'^ia^), ro^ro(0, 7o ^ 7o(fT'^) ■ (3.6) 

By slowly varying we mean that the derivatives of the world-volume fields are sufficiently small. 
In order to quantify this, we introduce a set of re-scaled coordinates cr" = eo"", e <^ 1, and 
consider the w.v. fields to be functions of a^. In this way each derivative will produce a factor 
of £. Moreover, two derivatives will be suppressed by a factor of e compared to one derivative 
and so on. Effectively what we are doing is to consider arbitrary varying world-volume fields (no 
restrictions on the size of derivatives) and "stretching" them by a factor of 1/e ^ 1. In this 
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way we will only consider slowly varying fields and the derivative expansion is controlled by the 
parameter e.^ The fields can now be expanded around a given point V 



V 



(3.7) 



We now seek derivative corrections to the metric and gauge field denoted by respectively ds'^^^ 
and ^(1), so that 



dsL + edsl-. + ©(e^) and A = A(o) + e^(i) + ©(e^ 



(3.8) 



solves the equations of motion to order e. By a suitable choice of coordinates, we can take the 



point V to lie at the origin a'^ = (0, 0). Moreover, we can choose coordinates so that u'" 



1(0,0) 



1, 



u 



(0,0) 



0, i = 1, ...,p (the rest frame of the boost in the origin) . In these particular coordinates, 



the O*'^ order metric ds^Q^ takes the form 



-2h 2 d-ydr 



4 = 1 



h r" 
+ 2 



+ 2iB 



1 ) a"'daUidvda 
1 



a^daUi do-Mr 



N 

2 j /i^/2+i r 



+ 



h \ rQ 



'(n+l) 



(3.9) 



where we have denoted ro|(o,o) = '"o and 7o|(o,o) = 7o- Clearly the system has a large amount 
of gauge freedom. Following the discussion of the definition of r^,, we want the r coordinate to 
maintain its geometrical interpretation. We therefore choose 



0, 



(3.10) 



and we moreover take 



9hn 



and ^1^) = . 



(3.11) 



The background g^g) exhibits a residual SO(p) invariance. We can use this to split the system 
up into sectors of SO(p). The scalar sector contains 4 scalars, a!^\ g^) , Qm and Tig^^^ . The 
vector sector contains 3 vectors fffjP and g^]^ . Finally, the tensor sector contains 1 tensor 



9^^ 



^ij^ — ^Tiglj^ (the traceless part of 



^In the end of the computation, we of course set e — 1 and keep in mind that tlie expressions only hold as a 
derivative expansion i.e. for sufficiently slowly varying configurations. 
*In these coordinates u" = 1 + 0{£^). 
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We parameterize the three SO{p) sectors according to 



Scalar: = -y^Nj^il + joY^h-'a,, gg) = /i^-^/^.r, gi'J = /I'V™, Tr^^ = /i^Tr/,,- , 

Vector: A« = -^iV^o(l + 7o) a., 9i]^ = h''U, 5^ = /i^-^/2/w , (3.12) 

Tensor: = h^'Jij , 

where f^j = fij — |Tr/jj. The parameterization is chosen in such a way that the resulting EOMs 
only contain derivatives of fab and Ua and will thus be directly integrable. 

3.2 A digression: Reduction of Einstein-Maxwell theory on 5""+^ x W~^^'^ 

Here we explain how it is possible to treat general n and p by integrating out the transverse 
non- fluid dynamic directions. 

In order to work out the full set of solutions and find the general form of the stress tensor 
and current, it is enough to consider fluid dynamic fluctuations in 1 + d (2 < d < p) directions of 
the brane.^ In particular, it is enough to consider d = 2. Indeed, since the background is SO(p) 
invariant, the correction ds^^^^ will consist of SO(j>) invariant tensor structures. The same holds 
for the effective blackfold stress tensor and current. In order to identify these tensor structures, 
it is enough to consider fluctuations in only 1 + d directions (time + d flat spatial directions) of 
the brane. Considering only fluctuations in 1 + d brane dimensions, the metric is of the form 
(reduction of the p-brane with n + 2 transverse dimensions) 

i=d+l 

with the one- form gauge field of the form Af^ = Aa{crf). Here the subscript / means 'fluid' since 
the d + 2-dimensional base space with the metric ds^^^ will contain the fluid dynamical degrees of 

freedom in our computations. Integrating out the S"'^^ and M*'"'^"^ (see appendix A), the EOMs 
of the system take the form 



R'£ = Tah + (n + 1) (VaV'VbV' + VaVfeV) + (p " C^) (V^^Vft^ + V^Vfe 

□V' + \(j> - d)Vb(t) + (n + l)Vb^] VV = ™e"2^ + K , 
U(t> + [{p - d)Vb(p + (n + 1)V6^] VV = « , 

VaF'^' = / , 

where the tensor J^ab^ vector and scalar k are given by 



(3.14) 



1 TP Tpoh 

T\ = -F'^'Fbc - k6%, f = F'^" ((n + 1)V,V + (p - d)V,0)) , k = °^ ■ (3.15) 

z 4(p + n + 1) 



Working with these effective EOMs allows us to treat a general number of transverse and brane 
dimensions. 



^We thank Joan Camps for pointing this out. 
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4 First order equations 



In order to compute the effective stress tensor and current and thereby extract the transport 
coefficients, we need the large r asymptotics of the perturbation functions which are decomposed 
and parametrized according to equation (3.12). We denote the first order Einstein and Maxweh 
equations by 

R,. - \f,,F/ + = e£,. + 0{e-) = , ^^^^^ 

VpFP^, = eM^ + 0{e^) = {). 

In this section we will find the solution to each SO(p) sector in turn and explain how the regularity 
on the horizon is ensured. 



4.1 Scalars of 50(jo) 

The scalar sector consists of seven independent equations which correspond to the vanishing of 
the components: £yy,£rv,£rr,'^'^£ii^£Q.Q,-,-M.v,-M.r-^ 



Constraint equations: There are two constraint equations; 8"^^ = and Mr = 0. The two 
equations are solved consistently by 

(n + 1 +p57o) 9„ro = -ro(l - B-^o)^^u' , (4.2) 

and 

(n + 1 + pBjo) 9„7o = -270(1 + -fo)diu' . (4.3) 

The first equation corresponds to conservation of energy while the second equation can be inter- 
preted as current conservation. These are equivalent to the scalar conservation equations given 
by (2.14) in the rest frame. 

We now proceed to solve for the first order correction to the scalar part of the metric and gauge 
field under the assumption that the fluid configuration satisfy the above constraints. Imposing 
the constraint equations will make £rv and 8^^ linear related and one is therefore left with five 
equations with four unknowns. 



Dynamical equations: The coupled system constituted by the dynamical equations is quite 
intractable. One approach to obtaining the solution to the system is to decouple the trace function 
Tr/ij- Once Tr/jj is known, it turns out, as will be presented below, all the other functions can 
be obtained while ensuring that they are regular on the horizon. 

It is possible to obtain a 3rd order ODE for Tr/jj by decoupling it through a number of 
steps. One way is to use £„■ to eliminate fl^ and then take linear combinations of the remaining 
equations. The resulting combinations can then be used to eliminate /^^ and f'^^ such that one is 
left with two equations in terms of and Tr/jj which can then be decoupled by standard means. 
The resulting equation is schematically of the form 

H^i^'^Xr) [T^U^f (r) + H^^^^\r) [Tr/,,]" (r) + <'^)(r) [TY/,,]' (r) = S^Ar) , (4.4) 

where Hi,H2 and H3 do not depend on the sources (world- volume derivatives) and the source 
term S'xr only depends on the scalar diU^. The expressions for these functions are however very 
long and have therefore been omitted. After some work, one finds that the equation is solved by 

Tr/,,(r) = 4^ + 7oc?,^G(r) - 2idiu')Trflf (r) , (4.5) 

''In the reduction scheme outhned in section 3.2 we have Sna — where £^ is the EOM for ^ given in (3.14). 
Similarly, we have TrSij = Trf"^ - (p - d)h'^£4,, where £^ is the EOM for 0. 
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where the terms containing the two integration constants c^^ and d^^ correspond to the homo- 
geneous solution. The entire family of homogeneous solutions to equation (4.4) of course has an 
additional one-parameter freedom which has been absorbed in the particular solution TT:f-j\r) 
and been used to ensure horizon regularity'^. The function G is given by 



The particular solution which is regular on the horizon is given by 



Tr4'^(0 = -^(1 + 7o)^«7oG(r) + ( r, - ^ (1 + 70)^ log /(r) ) (1 - P^oG{r)) 



n 



with the coefficients 

a = 2B 



2{n + l) +pB-fo 
{n + l+ pB-f^y 



and 13 = B 



n + 2+ pBjQ 
n + l+ pBjo 



With Tr/jj given, the equation £rr = will provide the derivative of frv, 



d 



2(n + l)+p5^7oU(r 



dr 



Mr)-[Tr/,,]'(r) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



Since the equation is a 1st order ODE, the regularity of the horizon is ensured by Tr/y. Note that 
it is possible to perform integration by parts and use that the derivative of takes a simpler form. 
One can thereafter obtain an analytical expression for the resulting integral. This expression is 
however rather long and does not add much to the question we are addressing for which we are 
in principle only interested in the large r behavior given by 



frvir) 



„(2) G{r)h{r) 

7oCTr 7 ;;r— 

2+pS^7o 



°o nk 

I I 



(4.10) 



k=l 



The first two terms constitute the homogeneous solution and the particular solution is given in 
terms of the coefficients c 
are given in appendix B. 



terms of the coefficients ai^ and /3™^ which depend on n,p, and 70. The first set of coefficients 



Using the expression for /^^ in terms of Tr/jj, the Maxwell equation = becomes a 2nd 
order ODE for the gauge field perturbation. 



d_ 

dr 



1 



nr 



d 



2(n + l)+pB^7o 



dr 



1 



[Tr/.,]'(r) 



(4.11) 



This equation is solved by a double integration. The inner integral is manifestly regular at the 
horizon, one can therefore work directly with the asymptotic behavior of the right-hand side 
before performing the integrations. The large r behavior of the perturbation function is thus 
found to be 



av{r) 



(1) « , (2) 



1 



c>, 'r 



+ 4'^ + ^7oc£^G(r) + (9, 



n 



I r 



n 



-r+y^ 



k=\ 



(4.12) 



''Note that equation (4.4) has been derived under the assumption that diV^ 7^ 0. This especially means that 
when there are no sources the one-parameter freedom disappears in accordance with (4.5). 
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where the first three terms constitute the homogeneous solution and the particular solution is 
given in terms of the coefficients al*^^ and /J^,'^^ depending on n,p, and 70. The first set of 
coefficients are given in appendix B. 



The last perturbation function fy^ can be obtained from 
f'j.y in terms of Tr/jj the equation is schematically of the form 



0. Using the expression for 



d 
dr 



= Gi [Tr/,,(r)] + G2 K(r)] + S,,(r) 



(4.13) 



where Gi, G2 are non-trivial differential operators and the source Su depends on diU^. Again, the 
full expressions have been omitted and we only provide the large r behavior, 



' r 



k=l 



(4.14) 



with the homogeneous part given by 



.(1) 



(2) „ 



2(l + 7o)-|(ci^)-cWro"7o) 



G(r) 



l + 27o--^7o|cS) 



(4.15) 

The solution is ensured to be regular at the horizon. The coefficients and pit^ depend on 
n,p, and 70. The first set of coefficients are listed in appendix B. 

Finally, one must ensure that the equations coming from TrSa and the angular directions 
i^nn = 0) are satisfied. This will impose the following relations between the integration constants, 



.(1) 

.(2) 



((n + p)cg) + 8(1 + 7o)(4'^ - 4iVo"7o)) • 



(4.16) 
(4.17) 



This completes the analysis of the scalar sector. The remaining undetermined integration con- 
stants are thus: , c^^ , Crv, oh'^ , for which Crv and cij^ will be fixed by requiring the spacetime 
to be asymptotically flat while the rest constitute the freedom of the homogeneous solution. Note 
that the above functions reproduce the neutral case as 70 (cr") — )■ 0. 



4.2 Vectors of SO{p) 

The vector sector consists of 3p independent equations which correspond to the vanishing of the 
components: £ri-,£vi and M.i. 



Constraint equations: The constraint equations are given by the Einstein equations 8"^^ = 
and are solved by 

diro = ro(l + iV7o)5^,Ui , (4.18) 

which are equivalent to conservation of stress-momentum. These are part of the conservation 
equations given by (2.14) in the rest frame. Similar to above we now proceed solving for the first 
order corrections to the metric and gauge field under the assumption that the fluid profile satisfy 
the above constraint (4.18). 
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Dynamical equations: The remaining equations consist of p pairs consisting of one Einstein 
equation £yi = and one Maxwell equation Mi = 0. The structure of these equations is the same 
as in the scalar sector. The Einstein equation £yi = is schematically of the form, 



while the Maxwell equation Mi = is, 

Mf^)(r)<(r) + Mt^'hMir) + M^^)(r)/^(r) = 5,(r) . 



(4.19) 



(4.20) 



Again the functions and M/., k = 1, . . . , 3 have been omitted. 

To decouple the system we differentiate £yi once and eliminate all aj(r) terms in Mi. Doing 
so, one obtains a 3rd order ODE for fvi{r) which can be written on the form. 



dr 



with 



1-Cl^ 



Cl 



d 
Ik- 



^ —f'.(r) 



V 



N -I 



-70 • 



(4.21) 



(4.22) 



l + iV7o 

It is possible to perform the first two integrations analytically and ensure regularity at the horizon. 
The first integration is straightforward while the second involves several non-trivial functions. The 
large r behavior of the fm function is found to be 



fviir) 



1 



fir) 
h{r)^ 



I n (fc) 



fc=l 



(4.23) 



where the first two terms constitute the homogeneous solution. The first set of coefficients a 

(k) 

and f3^,- are given in appendix B. Notice that the sum vanishes in the neutral limit. 
Once the solution of is given we can use to determine aj, 

+ 



(fc) 



a- 'r + j3l Wq 



(4.24) 



k=l 



where the first two terms correspond to the homogeneous solution. The first set of coefficients 

(k) (k) 

a- and /3f' are given in appendix B. 

The remaining undetermined integration constants are thus: c^^\c^i\ and c^^ . The constant 

.(2) 



c^- corresponds to an infinitesimal shift in the boost velocities along the spatial directions of the 

brane while c^^^ is equivalent to an infinitesimal gauge transformation. The last constant c^^^ will 
be determined by imposing asymptotically fiatness at infinity. 



4.3 Tensors of SO{p) 

There are no constraint equations in the tensor sector and p{p + l)/2 — 1 dynamical equations 
given by 

£ij - ^TriSij) = . (4.25) 



This gives an equation for each component of the traceless symmetric perturbation functions fij, 



d 



^K+V(r)4,(r)]=-a, 



'^0 



:r-(2(n+l)+pS-^7o)/i(r)^ 



(4.26) 
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where ^ 

aij = d(^Uj) - -Sijdkv!" . (4.27) 

The solution is given by, 

h{r) = Cij - 2aij {r, - l2 (1 + 70)^ log /(r)) , (4.28) 

where horizon regularity has been imposed and the constant Qj is symmetric and traceless and 
will be determined by imposing asymptotically flatness. 



4.4 Comment on the homogeneous solution 

We have now obtained the solution to the Einstein-Maxwell equations for any first order fluid 
profile which fulfill the constraint equations. These have been provided in large r expansions 
and are ensured to have the right behavior at the horizon for any of the remaining integration 
constants. One remark that is worth mentioning is that fri did not appear in the analysis above 
and corresponds to a gauge freedom. This gauge freedom does not play a role for n > 2, but is 
expected to play a role for n = 1 to ensure asymptotically flatness. 

We now want to provide some insight into the meaning of the remaining integration constants. 
One can separate the constants into two categories; the subset that are fixed by asymptotically 
flatness and the subset that corresponds to the e-freedom of the parameters in the zeroth order 
fields. The latter corresponds exactly to the remaining freedom of the homogeneous solution. In 
the above the homogeneous part of the fields are given exact. 

One finds that the homogeneous part of the scalar sector corresponds to shifts in tq — )■ ro+e5ro, 
7o — ^ 7o+£<^70) and the gauge freedom — )• a^+eday of the zeroth order metric given by equation 
(2.6). Indeed, performing the above shifts and redefining the r coordinate, 

r r(l- c- J^ ri5logro + 6logjo \ 
\ r"- n + ph{r) J ' 

such that the angular directions does not receive first order contributions in accordance with 
the gauge choice (3.11), one can relate the integration constants to the two shifts and gauge 
transformation by, 

c^^ = 2B{n6 log ro + 6 log 70 ) , 

c(,2) = -n6 log ro - J; +^^0 a log 7^ ^° Sa^ ^ (4.3O) 

2(1 + 70) \/A^7o(l + 7o) 



5ay 



r^V^A^7o(l + 7o) 



For the vector sector one finds that the homogeneous part corresponds to the shift of Uj — ?• 
Ui + edui and the gauge transformation — )• + edai . The first transformation corresponds to 
global shifts in the boost velocities. In the same r-coordinate, one has 

cf^ = - , ,. (4.31) 

ViV7o(l + 7o) 

This accounts for all the e-freedom in the full solution. 
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4.5 Imposing asymptotically flatness 



We now turn to imposing the boundary condition at infinity, namely requiring the solution to be 
asymptotically flat. To impose this we must first change coordinates back to the Schwarzschild- 
like form. Moreover, we need the fields expressed in Schwarzschild coordinates for obtaining the 
effective stress tensor and current. In order to change coordinates, we use the inverse transfor- 
mation of the one stated in equation (3.1). The transformation can be worked out iteratively 
order by order. To first order the transformation from EF-like coordinates to Schwarzschild- like 
coordinates for a general rQ^a"") and 70(0"") is given by. 



t + r^ + e 



(t + r^) {dror^dtro + d^gr^dao) + {drarAro + d^^^rAlo) 



a 



x' + e 



u 



(4.32) 



r. + 0(e'). 



It is now possible to transform all the fields to Schwarzschild coordinates and impose asymptoti- 
cally fiatness. This leads to 



0, 



0, 



-Tr 



0, 



0. 



(4.33) 



We now have the complete first order solution for the black brane metric and Maxwell gauge field 
that solves the Einstein-Maxwell equations. 



5 Viscous stress tensor and current 

In this section we will compute the effective stress tensor and current of the first order solution 
obtained above. Before doing this, we shall briefly discuss the general form of the first order 
derivative corrections to the stress tensor and current. 

5.1 First order fluid dynamics 

We write the stress tensor and the current as 

T<^' = T^^^+Ilfl^ + 0{d^), f=j-^^ + Tl^^+0{d^), (5.1) 

where the perfect fluid terms were written down for our specific fiuid in section 2.2. The tensors 
n"-|'^ and T^j'j are the first order dissipative derivative corrections to the perfect fluid stress 

tensor and current, respectively. The speciflc form of 11"^^ and T^j*^ are encoded in the first order 
correcting solution obtained in the previous section. As is well-known, to any order in derivatives, 
it is in principle possible to write down all the terms that can contribute to the stress tensor and 
current (see e.g. [32]). In this way the dissipative corrections to the stress tensor and the current 
can be characterized in terms of a set of transport coefficients. It is possible to show that the 
most general form of 11"]'^ is given by® 

n^^-^ = -2i]a''^ - C^A"'' , (5.2) 

where cr'^'' is the usual shear tensor and is given by 

= A'^^ - Ae/) A'^" with ^ = dau\ (5.3) 

*Here we have imposed the Landau frame gauge on the stress tensor Uall"*) = 0. Similarly one can impose a 
Landau frame condition on the current. It takes the form «.aT?j) — 0. 
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The coefficients ij and C are respectively tlie sliear and bulk viscosity transport coefficients and 
were computed for the neutral brane in [4]. The bulk and shear viscosities are associated with 
the scalar and tensor fluctuations, respectively. Note that although the overall form of n^-j*^ is 
the same as in the neutral case, the transport coefficients are now expected to depend on both 
the temperature and the charge i.e. on both tq and 70. Also note that the viscosities r/ and ( are 
required to be positive in order to ensure entropy creation in the fluid [33] . 

Using similar reasoning, it is possible to show that the most general form of T"-^^ (in the 
Landau frame) is given by^ 

^W = -^(^)'^'''^^(|) • (5.4) 

Here S is the charge diffusion constant which is associated with the vector fluctuations. Indeed, 
it is possible to derive that with D > 0, the term (5.4) is the only term which can be constructed 
from the fields and that is consistent with the 2nd law of thermodynamics [33]. Plugging in the 
values of ^ and T in terms of rg and 70 and using the vector constraint equation (4.18), we find 
that (in the rest frame) 

T^i)=0, Tji)~7o(l + 7o)5.^x' + ^5ao. (5.5) 

Since the derivatives appear in a very specific combination in this expression, this in fact provides 
us with a non-trivial check of the blackfold fiuid description. 



5.2 Computing the effective stress tensor and current 

The quasi-local stress tensor r^^ is obtained by background subtraction. We consider a surface 
at large r (spatial infinity) with induced metric h^i, and inwards pointing normal vector and 
compute the components of the quasi-local tensor by, 



SvrGr, 



flU 



flU 



hfMuK 



K, 



fJ,U 



(5.6) 



where K^^, is the extrinsic curvature of the surface and K = g^^K^^. The hatted quantities are 
the subtracted terms which are computed on flat spacetime with the same intrinsic geometry on 
both boundaries. Notice that the transverse space bear the structure dQ'^^_^-^^y One finds that 
for the transverse directions r^n = while for the brane directions we obtain the fluid stress 
tensor by, 



hm^ 

r— >oo 2 



r Tab- 



(5.7) 



where is the volume of the (n + l)-sphere. We find 



~ ~ levrG 
+ ero"+i(l + 7o)^ 



Sij (ro + e(<5ro + x^daVo))" 



2 a, 



P 



5ijdku'' 



2 (n+p+l)(n + l) , 
H ; — To- "ij Ok u 



(5.^ 



®It is possible to include a parity violating term as was found in [23]. However, since we have no Chern-Simons 
term in the theory such a term is not relevant. 
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where the expressions are vahd to order 0(e). In a similar manner the current is obtained from 
large r asymptotics of the gauge fields. Ensuring that the Lorenz gauge condition V^A^ = is 
satisfied, the current is obtained using 

J, = hm "^^^r-A, . (5.9) 

One finds 

Jt = J^^^ (^0 + e{5ro + x«5aro))" v^Toll + 7o) + ei^lo + x-da7o){l + 270) ' , 

, ^(n+l) rj-p n HTT^ v( fx , \ lojl + 7o)^^,|3^ + ^dgp \ ^^'^^^ 

-^i = ' / WA^ro V7o(l + 70 e{5uj + x daUj) - ero ^ . 

16^^ V n(l + iV^o)7o(l + 7o)^+V 

Again these expressions are valid to 0(e). It is now possible to read off the transport coefficients. 
Before doing this, we require that the Landau frame renormalization conditions 11**^ = 11*^^ = 
and T*-|^^ = are satisfied. Equivalently we require the shifts 6rQ and 6jo of the zeroth order 
solution to vanish. Notice that the stress tensor and current do not depend on the gauge trans- 
formation 6aa as they should of course not do. Also recall that the shifts were related to the 
integrations constants by (4.30). 

Setting 5rQ = ^70 = 0, the shear and bulk viscosities are determined using the form given by 
equation (5.3), 

„_%^,nH-i^ , C 2 (n + p + l)(n + l) 

The second term of Ji is seen to have the right proportionality according to (5.5) and hence using 
the form of equation (5.4) the diffusion constant can be determined, 

® _ 1 _ 

Notice that all the transport coefficients are found to be positive which is expected for a consistent 
effective fluid dynamic theory. We have now obtained the first order derivative corrections to the 
effective stress tensor and current. 

5.3 Hydrodynamic bounds 

We will now check the result of the shear viscosity against the expectation that the transport 
coefficient should satisfy the bound 

2>J-. (5.13) 

s 47r 

Using equations (2.9) and (5.11), the system is seen to saturate the bound as expected. 

In addition, it is worth to investigate the bulk to shear viscosity ratio proposed by ref. [16], 

C^^/'i ,2 



where Cg is the speed of sound. Although one should keep in mind that the proposal of this bound 
relies heavily on holographic considerations, we find when using the value given by equation (6.5) 
for the Reissner-Nordstrom brane, that the bound is satisfied in the range 



„ , . n+1- ^l + n(n+p + 2)' 

< 70 < , (5.15) 
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while for large values of 70 the bound is found to be violated. If we instead of Cs in (5.14) use 
the proposed quantity [17, 18] 



(f), 



1 



1 + n 



1 + 270 
1 I PB 



(5.16) 



computed for fixed charge density Q, we find that the bound will always be violated (except for 
the neutral case where cq = Cg). 



6 Stability and dispersion relations 

In ref. [3] the Gregory-Laflamme instability was successfully identified with the unstable sound 
mode of the neutral black brane. This analysis was further refined in [4] and considered for branes 
charged under top-form gauge fields in [31]. In this section we address the issue of stability and 
dispersion of long wavelength perturbations of the Reissner-Nordstrom black brane. Moreover, 
we comment on the connection to thermodynamic (in)stability. 

6.1 Dispersion relations 

It is straightforward to show that the first order fluid (conservation) equations take the form 

e = -{w - C'&j'd - 2r]aabcr , u = — ^ ^ 

w — Cv , , 

2 N (6-1) 

where the transport coefficients and the factor associated to 3 are coefficients in the derivative 
expansion and should be treated as constants. In order to find the speed of sound and dispersion 
relations, we consider small long wavelength perturbations of the fluid 

$^$ + (^$e<'^*+'=^^'), r^r+(5re<'^*+^'^^'), = (1,0,...) ^ (l,(5'a^e<'^*+'=^^')). (6.2) 

The charge density Q, energy density e, and pressure P are perturbed according to 

Q^Q + ,^Qe^('^*+^'^^'), e ^ e + 5ee^H+'=^^'), P^P + <^p e<'^*+^'^^') , (6.3) 

where the amplitudes can be expressed in terms of thermodynamic derivatives that depend on 
the specific equation of state. Note that 5p = Q6^ + s6T as a consequence of the Gibbs-Duhem 
relation. Plugging the expressions into the first order fluid equations (6.1) and linearizing in the 
amplitudes, we obtain the p + 2 equations 



( ) + (If ) + + (1*^ - ») = ° ■ 

We stress that the thermodynamic derivatives are not dynamical and do only depend on the 
equation of state of the fluid in question. In our case they can be computed from (2.9) and (2.12). 
In order to find the w that solves this system for a given wave vector /c* , we set the determinant of 
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70 



Figure 1: The qualitative behavior of the sound mode uj = Cgk + ak"^ + 0{k^) given by equation 
(6.8) as a function of 70. The leading order correction (speed of sound) is seen to be positive 
when the charge density passes the threshold 70 = (-D — 3)/2 (indicated with the vertical dashed 
line) while the first order correction diverges at the threshold and becomes negative. The two 
horizontal dashed lines correspond to the speed of sound squared and first order correction in the 
neutral case, respectively. 



the system of linear equations in the amplitudes to zero. To linear order in /c* (i.e. at the perfect 
fluid level) the dispersion relation gives the speed of sound Cs = uj/k. Using the equation of state 
(2.12) and solving the system to linear order, one finds 

d^i'-f] =-i:i|2»(„+l+pB,„)-\ (6.5) 



9e y ^ 1 + iV7o 

As was found with the p = q branes of supergravity [31], the speed of sound only depends on 
the charge parameter 70. For zero charge 70 = we recover the neutral result = — l/(n + 1). 
Since a negative speed of sound squared signifies an unstable sound mode, the neutral brane is 
unstable under long wavelength perturbations. Indeed, this instability is exactly identified with 
the GL instability [4]. However, as we increase 70 the speed of sound squared becomes less and 
less negative and for 

70 > ^ , (6.6) 

the q = brane becomes stable under long wavelength perturbations to leading order. Notice that 
the condition (6.6) can be satisfied for any non-zero charge density if the black brane temperature 
is low enough. Indeed, stability is obtained for T~ (GQ)~^/" (where the exact numerical factor 
depends on the number of transverse and brane dimensions). 

In order to check stability to next to leading order, we now work out the dispersion relation 
for the fluid to quadratic order in k. We solve the system of equations to 0{k'^). Solving for the 
longitudinal modes, we find the equation 



^_cl--t-{2{l--\r^ + Q\+—^(ni{-\\^]+0{k^) = Q, (6.7) 
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where the coefficients TZi, TZ2, and TZ (introduced below) are given in appendix C. Solving for the 
sound mode(s), we find the dispersion relation 



w(fc) = ±csk + — 1^ - - j + I j + ik^WZ . (6.8) 

For a general fluid both the first order term (c^) and the second order term must be positive in 
order for it to be dynamically stable. In this case, the above equation describes dampening of 
the long wavelength sound waves in the fluid. However, for our specific fluid, this can never be 
attained. Fig. 1 shows the general behavior of Cg and the (second order) attenuation term in 
(6.8). The conditions on 70 for dynamical stability are found to be complementary; when the 
sound mode is stable to leading order it is unstable next-to-leading order and vice versa. The 
Reissner-Nordstrom brane thus seems to suffer from a GL instability for all values of 70. 
In addition to the sound mode we have a longitudinal diffusion mode given by 

^{k) = —±k^ = -—— —— . (6.9) 

^ ' clw 47rr(l - S70) ^ ^ 

We see that in general this mode is stable if and only if T^-i/c^ > 0. In our case this amounts 
to the condition 70 > -^^^^ which is the same as (6.6). Finally, we also have a shear mode which 
takes the form 

oo{k) = '-^k\ (6.10) 

w 

The fluctuations of the shear mode are very simple, they are transverse displacement of effective 
fluid with no variations in the charge and energy densities. Notice that this mode is always stable. 

6.2 Thermodynamic stability 

The conditions for thermodynamic stability of the Reissner-Nordstrom black brane are computed 
in the grand canonical ensemble since charge is allowed to redistribute itself in the directions of 
the brane. Using the thermodynamic quantities in equation (2.9) and (2.12), one finds the specific 
heat capacity C and the (inverse) isothermal permittivity c to be, 

C=lL'\ _ + 1 + {2- n{N -2))jo 



TJq V (niV-2)7o-l , 
$\ / 1 \ 1 ^'-l') 



QJt V(7o + l)(l- (niV- 2)70) y sT' 

Thermodynamical stability is obtained if the two quantities are positive. However, these two 
conditions are complementary and can never be satisfied. This is also what was found for the class 
of smeared Dp-branes considered in e.g. [21]. Indeed, this complementary behavior is analogous to 
what was found for the dynamical analysis. However, the critical value of 70 where the quantities 
switch sign is not coinciding for the two analyses. It would be interesting to further investigate 
how the instability predicted by the dynamic analysis and the thermodynamic computation are 
related thus making a more precise connection to the correlated stability conjecture in the charged 
case [21]. 
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7 Discussion 



We have investigated the nature of the hydrodynamic effective theory that governs the intrinsic 
long wavelength fluctuations of the Reissner-Nordstrom black brane. Our analysis has extended 
the established cases of the interrelation between gravity and fluid dynamics. Although the 
analysis of section 4 is quite technical, the problem at hand provides the purest example of a 
black brane carrying charge. With the extraction of the effective stress tensor and current, our 
analysis has provided the generalizations of the known neutral shear and bulk viscosities. We 
find that the shear viscosity receives the expected modification such that i]/ s = 1/Att. Note that 
the entropy has the form as given in equation (2.9) for the entire family of generalized Gibbons- 
Maeda black branes, we therefore expect the result for r/ given by (5.11) to hold in general. In 
particular, this includes the case of the D3 brane. The bulk viscosity was found to be non-zero 
positive for all values of the charge as expected since the effective fluid is not conformal. The Q/ri 
bound proposed by ref. [16] was found to be violated for certain values of the charge parameter 
while it was demonstrated to violate the bound proposed in [17] in the entire range of non-zero 
7o, thus providing a counter-example. Finally, we computed the charge diffusion constant of 
the Reissner-Nordstrom black brane. We note that, as with the shear viscosity r], the value of 3 
given in (5.12) only depends on N which could be an indication that the result will hold for more 
general cases where e.g. the black brane is charged under higher form gauge fields. 

The speed of sound was found to be imaginary for small charge densities, but becomes real for 
sufficiently large charge parameter 70 > {D — 3)/2. For large charge density it therefore seems that 
the Reissner-Nordstrom black brane is GL stable under long wavelength perturbations. However, 
including the first order corrections to the dispersion relations, one finds that the attenuation of 
the sound mode becomes negative above the threshold value of 70. The Reissner-Nordstrom black 
brane is therefore GL unstable for all charge densities, although it is worth noting that the brane 
is "less" unstable above the threshold. This complementary behavior of the instability is also 
reflected in the thermodynamic stability analysis where the specific heat capacity and isothermal 
permittivity show a similar behavior. It would be interesting to investigate the relation between 
the two approaches in more detail, that is, establish a more precise connection to the correlated 
stability conjecture [21]. Also, it would be interesting for comparison to perform a numerical 
analysis of the long wavelength perturbations in the current setting as was done in the case of 
the neutral brane, where excellent agreement was found. 

Regarding the stability analysis, it is also worth noting that since the value obtained for the 
bulk viscosity leads to violation of the C/rj bound, one might question the validity of the stability 
analysis for the case of a black brane charged under a top-form gauge field examined in ref. 
[31]. Here the dispersion relations were written down using the assumption that the C/rj bound 
proposed by [16] is saturated. 

An interesting computation, that has not been investigated in the blackfold literature, is the 
computation of the entropy current a la [34]. Computing the entropy current could provide 
a consistency check of the transport coefficients and the framework. We hope to address this 
question in the future. 

Another natural future direction of this work is the generalization to black branes charged 
under higher form gauge fields possibly in the presence of a dilaton field. Of particular interest to 
string theory this would include the black Dp-branes that carry charge under a Ramond-Ramond 
field for which the case of the D3-brane would be included. The case of the D3-brane would also 
be interesting in the context of the AdS/CFT correspondence. This could namely help elucidate 
possible relations between the AdS and flat space case i.e. the connection between the blackfold 
approach and the fluid dynamical regime of AdS/CFT. In this regard, it would also be interesting 
to understand, in a systematic manner, how the map of ref. [28] extends when matter fields are 
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included. This could also work as a method for obtaining the second order transport coefficients. 
However, taking the computation to second order from first principles (as in this work) is also of 
interest. 

Finally, it would be interesting to include a Chern-Simon term in the theory. This was 
considered in AdS fluid/gravity in the papers [22, 23, 35]. However, we note that black brane 
solutions analogous to the generalized Gibbon-Maeda solution with such a term in the action is 
to our knowledge not known in the literature. 
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A Reduction 



In the first part of this appendix we will show how the equation of motions for the general case 
of a reduction of an Einstein-Maxwell theory on an Einstein manifold can be obtained. In the 
second part we will provide the example of applying the procedure for d = 2 on the zeroth order 
solution. 



A.l Reduction of Einstein-Maxwell theory on an Einstein manifold 

We consider Einstein-Maxwell theory on a D-dimensional space of the form 

ds2 = g^Ax'^dx" = dsl^ + e^'^^^'-^ds^^;) . (A.l) 

Here ds^^^ denotes the metric of the base manifold M.{h), x*^^^ denotes the coordinates on M{b)^ 
■0 is a function on M{b) and ds^^^ is the metric of an Einstein manifold Mf^E) with coordinates 
x^^y Since M.(^e) is an Einstein manifold, we have 

dER^''^ = Re9^''\ (A.2) 

where dg, g^^\ R^^^ and Re are respectively the dimension, the metric, the Ricci tensor and 
(constant) curvature scalar of M[e)- Moreover we consider a gauge field (minimally coupled to 
gravity) which only depends on x^^f^-^ and only takes values along the base manifold A4(^b)- 
Schematically 

A^{x) = Ai{xb) . (A.3) 



The action S of the system is given by 



S = Sg + SEM, Sg= I d^xy/\g\'R, Sem = / d^xy/\g 



(A.4) 



where R denotes the Ricci scalar of the full metric We can now perform a reduction and 
integrate out Ai(^E)j one finds 



Sg^ [ d^'^xtVHl'e''^^^'"'^ \Rb + REe^^^""^ + 4(V0)'| , 
5em~ / d'^^xtyWe^^^^'^'^FijF'^ . 



(A.5) 



Having worked out the reduced action, it is easy to work out the equations of motion. As usual, 
the resulting system is EM theory on A^f, coupled to a dynamical scalar field and a current. The 
EOMs are 

= 1 (if, ^Fju - + (VaV-VfeV + VaVbV) , 

V^F'^ = dEF^^Vill) . 
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A. 2 Reduction of the zeroth order solution 

In this section we demonstrate how the reduction works for the 0**^ order solution with (fluid) 
dynamics in two spatial directions (in other words, an ordinary boost in the ((Ti,(T2) direction). 
Now the base space is composed of the three fluid brane directions (one time and two spatial 
directions, ((7^,0"^) along with the radial direction r). The metric has the form 



(A.7) 



with a, b = 0,1,2 and where Xy , i = 3, ...,p are the p — 2 static brane directions. We now integrate 
out the transverse sphere and the p — 2 brane directions. The functions ip and (p a-re given by 



(f){r) = Tp{r) + 2 log r = B log h{r) . 



It is now straightforward to compute k, and J-"^. Here denotes coordinates of the four 
dimensional base space x'* = {a^,a^,a'^,r). One finds 

Bn"^ /ro\2n ^q(i + ^o) 



dx" 



n^'-Ub + 1 



/i^-i(r 
2 



N 



B n 



2('^)"7o)n'^5.„,, (A.9) 



' dr . 



It is now possible to show that, as expected, the reduced system obeys the EOMs with these 
effective sources. The above sources get derivative corrections in the perturbative expansion. 

B Coefficients of the large r expansions 

In this section, we list the first set of large r expansion coefficients of the metric and gauge field 
given in section 4. 



Scalar sector 

Below is listed the first set of coefficients for the large r expansions of f^.^. 



a, 



(1) 



n((n + pf + (n + p){2{p + 1) + n{p + 2))7o + 2p{p + 2)7^) 



(n - l)(n + p)2((n + 1) +PS70) 
Below is listed the first set of coefficients for the large r expansions of a^, 

(1) _ n{{2n + l)(ra + pf + (n + p){l + n{2n + 3)(p + l))7o + 2p(l -p + 2n{p + 1)7^) 



(B.l) 



a; 



/3i 



(1) 



(l + 7o)~ 
n 



1 - pjoB 



n - l)(2n - l)(n +p)2((n+ 1) +p57o) 
2(n + 1) +pBjo 



((n + l)+pB7o)2 



(B.2) 



Below is listed the first set of coefficients for the large r expansions of 

(1) 1 f n , n ^ , 47o(n+p(l+7o)) ^^ 

= + + (n + p)2((n+l)+pi?70))J 



{n 

/3W = -(l + 7o 



2(n + l) 



(B.3) 



n(n + 1 +pi?7o)^ 
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Vector sector 

Below is listed the first set of coefficients for the large r expansions of fy 



a 



a 



i-^ = Mho 

,(2) _ 



+ + + 1)(1 + 270) 



(n — l)(n + pY 



n + p + 1 
(n — l)(n + p) 



n+ p + 1 



2(n- l)(2n- l)(n+p)3 
2{n + p){n{n + p) + {An^ + n - 1 - 2p + 4np)7o)^ (5i7o) 

(^7o(4n(n + (n + - 2p + n(-3 + 2p + 4n(4 + 3n + 3p)))7o 

+ 4n(l + n + p){-l + n + 4n^ -2p + 4np)7^^ {d^f^i) 



(B.4) 



.(1) 



,(2) 



/ 270(1+ 70) (9^,^Xi) + (5,70) 



^™ V (l + 7o)-(l + A^7o) y 
Below is listed the first set of coefficients for the large r expansions of bi, 



a 



(1) 



2(n- 1) 



2{n + P + 2n{n +p + 1)70) 



n+ p 



(5„/3i) + 



1 + 270 
70(1 + 70) 



{dao) 



n + p 



(B.5) 



_(n + p + 1)70(1 + 70^ 

C Thermodynamic coefficients 



In this appendix we list a number of thermodynamic coefficients related to the analysis of section 
6. The two coefficients TZi and 7^2 are given by 



dQ 



dT 



Uo = -n^ 



T 



( de 



+ $ 



T 



id] 



Writing out the speed of sound given in equation (6.5) it takes the form 

7^l 



Q{Q 



Jde_ 



T 



Finally the coefficient associated to the dispersion relation of the sound mode is given by 



n 



1 7^? 



2 Q^w^c'i \ \dT 



Q 



f de_ 



T 



el-- 



dQ 



dQ 
dT 



T 



(C.2) 



(C.3) 



For the Reissner-Nordstrom solution we have 

7^l iV7o 7^2 l-A^7o(l + 27o) + n(l + iV7o)2 



r n + l+pS7o' sT^ 



1 + 270 + n(l - 570) 



n 



w 



2iV272(l + 7o)^ 



(C.4) 



sT^ {1 - Bjo){n + 1 + pB-fo) 
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